
Estimation of Surface Temperature and Heat Flux 
by Inverse Heat Transfer Methods using Internal 
Temperatures Measured while Radiantly Heating a 
Carbon/Carbon Specimen up to 1920°F 

Michelle Pizzo 

Intern, NASA Langley Research Center, Hampton, VA 
Applied Mathematics Ph.D. Student, Old Dominion University, Norfolk, VA 

Dr. Kamran Daryabeigi Dr. David Glass 

NASA Langley Research Center, Hampton, VA NASA Langley Research Center, Hampton, VA 


Math Awareness SIAM Conference 

Old Dominion University, Norfolk, VA 
April 11, 2015 


OUTLINE 


Problem Statement and Background 
Computational Code Development 

• Algorithm derivation 

• Data filtering 

• Mesh convergence 

• Method validation 

Analysis of LaRC Radiant Test Data 
Concluding Remarks 
Future Work 


PROBLEM STATEMENT 


• Purpose is to estimate surface temperature and heat flux values from internal 
temperatures measured while radiantly heating a carbon/carbon specimen up to 
1920°F 

• Initial and boundary conditions (BC) for T and/or q" must be known 


Ideally, BC's given on upper and lower 
surfaces of the material; direct 
problem needs only be solved 


upper BC, x = 0 

Material 

Direct 


lower BC. x = L 


In some cases, BC's given 
internally to material surface; 
both direct and inverse 
problems must be solved 


Material 


i 

upper BC, x = x 1 

i. 

Inverse 


Direct 

lower BC, x = x 2 



Lower TC 


Upper TC 


TC: thermocouple 
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COMPUTATIONAL METHOD 


DIRECT SOLUTION 


C p (x, T)p(x) 


dT(x, t) 

at 


a 

dx 


U(x, T) 


dT(_x, f)\ 
dx ) 


+ Q(x,t) 


• Developed a program to first solve the direct heat conduction problem 


- Direct Problem: solving between internal measurements 

• Created a mesh between the upper and lower data measurements 

• Calculated T and q" values at nodal points in the mesh and every time step using a 
one-dimensional, implicit, centered, finite volume method 
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COMPUTATIONAL METHOD 

INVERSE SOLUTION 


• Developed a program to then solve the inverse heat conduction problem 

- Inverse Problem: using direct problem solution to march to the surface of the material 

• Used work presented by A.S. Carasso, 1992 to select a space marching scheme 

• Space marching allows for the estimation of surface T and q" 

THE FUTURE 0 SPACE MARCHING TECHNIQUE WAS SELECTED FOR THE INVERSE ANALYSIS. 
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Carasso, A. S. "Space Marching Difference Schemes in the Nonlinear Inverse Heat Conduction Problem." 

Inverse Problems 8.1 (1992): 25-43. Print. 
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COMPUTATIONAL METHOD 

DATA FILTERING 



dT/dt vs. Time 

Noisy data compared with filtered data 



Time (s) 


dT(x, t) 

C v {x, T)p(x) = — [k(x, T) + Q{x, t) 

• In this study, local Blackman and Hamming windowed sine filters were used to 
reduce noise in the data. 

- Noise is reduced at each data point using information from a certain number of data 
points depending on the selected bandwidth 

• Bandwidth of 0.1 uses 41 data points, 20 on either side of each data point 

• Bandwidth of 0.2 uses 21 data points, 10 on either side of each data point 

- Saves on computational time when there is an excessive amount of data points 
compared to using a global filter 
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COMPUTATIONAL METHOD 

MESH CONVERGENCE 

• Solutions were calculated using multiple iterations of increasing nodal counts (3, 5, 
9, 17, 33, etc.) 
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(x — x 2 ) Iteration 1 Iteration 2 Iteration 3 


• Variance between T values at the current and previous nodal count was calculated 
over the entire range of time at the midpoint between the upper and lower 
boundary 

• Solution was considered to be converged once the variance was less than 0.01 K 


7 


METHOD VALIDATION 

COMPARISON TO EXACT ANALYTICAL SOLUTION 


• To solve for an exact, analytical solution the following assumptions were made: 

- Constant thermal properties dT(x,t ) d 2 T(x,t) k 

— ct cc — 

- No internal sources of heat dt dx 2 ' C p p 

• Spatial range: 0 < x < L 

• Constant initial condition: T(x, 0) = T 0 

• Constant boundary conditions: 7(0, t) = ^andTCL, t) = T 2 

• Solved for the unique analytical solution over a range of t = [0,60] seconds at 
varying spatial locations of x = L / 4 , L I 2 , and 3L / 4 
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r, W vuo i (m)_ * (^) p « £?] 


Test 1 

200 

100 

20 

0.01 

4 

1500 

1580 

1.69E - 

6 

Test II 
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Temperature Difference (K) Temperature (K) 


METHOD VALIDATION 

COMPARISON RESULTS 


Temperature vs. Time 
Test II, At = 0.05 seconds 



Upper boundary 

Converged: x = L/4 
4— Analytical: x = L/4 

Converged: x = L/2 
Analytical: x = L/2 
Converged: x = 3L/4 
Analytical: x = 3L/4 
Lower boundary 



1.0 

0.8 

0.6 

0.4 

0.2 



Temperature Difference vs. Time 
Test II, At = 0.05 seconds 


This temperature difference plot depicts 
the temperature difference between the 
analytical and computational solutions 
over the entire time domain. 
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METHOD VALIDATION 

COMPARISON RESULTS 



Heat Flux vs. Time 
Test II, At = 0.05 seconds 



Converged: x = L/4 
♦- —♦■Analytical: x = L/4 
Converged: x = L/2 
♦- —♦Analytical: x = L/2 

Converged: x = 3L/4 
♦- —♦Analytical: x = 3L/4 
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TC PLUG DESCRIPTIONS 


• Carbon/carbon (C/C) material, slightly different for each plug 

• For accuracy of solution 

- Material properties must be accurate 

- TC depths must be accurate 

- Analysis must represent the physics (ID analysis, ID physics) 


Horizontal-ply plug 


Vertical-ply plug 


SAME MATERIAL and 
SAME FIBER ORIENTATION 


DIFFERENT MATERIAL and/or 
DIFFERENT FIBER ORIENTATION 


LaRC RADIANT TEST SETUP 


• 6 in. x 6 in. x y 2 in. carbon/carbon test sample 
placed directly on Vi in. alumina insulation on 
top of water-cooled plate; heated with quartz 
lamps 

• Measured internal T using horizontal-ply and 
vertical-ply TC plugs, each with 4 embedded 
TC's; data sampled at 10 Hz 

• 12 different tests conducted with target T of 
500°F to 1920°F, Test 9 results discussed 


Front surface 
measurement 


Back surface 
measurement 


TC plug 





TCI 


TC2 


0.10 in. 
0.12 in. 
0.10 in. 
0.12 in. 


TC4 TC3 


Quartz lamp array 


V 4 in. Inconel septum plate 



V 8 in. C/C V 2 in. C/C V 2 in. flexible Water-cooled 
picture frame test sample alumina plate 
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Temperature (K) 


LaRC RADIANT TEST 9 

HORIZONTAL-PLY RESULTS 


TC plug 
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Front surface q" estimated by direct problem with front/back surface T as BC's 
Front surface T data, estimated q" assumed to be the correct surface values 
TC plug data analyzed with combinations of TCI, TC2, TC3, and TC4 as BC's 

- Surface estimations compared to the assumed correct values 

- Results indicated ACCEPTABLE inverse analysis T and q" surface estimations 



Temperature vs. Time 
Front Surface Comparison, LaRC Test 9 
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Heat Flux vs. Time 

Front Surface Comparison, LaRC Test 9 
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Error (%) 


LaRC RADIANT TEST 9 

HORIZONTAL-PLY RESULTS 


TC plug 


0.0% 

-0.5% 

- 1 . 0 % 

-1.5% 


Front surface q" estimated by direct problem with front/back surface T as BC's 
Front surface T data, estimated q" assumed to be the correct surface values 
TC plug data analyzed with combinations of TCI, TC2, TC3, and TC4 as BC's 

- Surface estimations compared to the assumed correct values 

- Results indicated ACCEPTABLE inverse analysis T and q" surface estimations 



Temperature Error % vs. Time 
Front Surface Comparison, LaRC Test 9 


Heat Flux Error % vs. Time 
Front Surface Comparison, LaRC Test 9 
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Heat Flux (W/cm 2 ) 


LaRC RADIANT TEST 9 

VERTICAL-PLY RESULTS 


TC plug 


• Data analyzed using thermal properties of the plug and of the 
surrounding material with TC2/TC4 as BC's 



TCI 

TC4 


TC2 

TC3 


Analysis using thermal properties of plug 
Heat Flux vs. Time 

Front Surface Comparison, LaRC Test 9 


Analysis using thermal properties of surrounding material 

Heat Flux vs. Time 

Front Surface Comparison, LaRC Test 9 



TC2/TC4 

front/back 


Time (s) 


Time (s) 


Analysis demonstrated a ID modeling technique can only consistently predict accurate surface 
T and q" when the TC plug properties and fiber orientation are the SAME as the surrounding 

material, hence approximating ID physics. 
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LaRC RADIANT TESTS 

CONCLUDING REMARKS 


• Developed computational code to solve ID direct and inverse heat conduction 
problems 

- Used one-dimensional, implicit, centered, finite volume method for direct problem 

- Used space marching techniques for inverse problem 

- Validated direct problem accuracy by comparing to an exact analytical solution 

- Accounted for data filtering and mesh convergence 

• Applied ID analysis to real experimental data 

• Demonstrated solution has better accuracy when using TC plugs manufactured 
with same material and fiber orientation (approximately ID physics) vs. when 
manufactured with different material and/or fiber orientation (non-ID physics) 

- Analysis yielded small errors for 7, q" when using horizontal-ply plug 

- Analysis yielded large errors for 7, q" when using vertical-ply plug 


Horizontal-ply 



Vertical-ply 
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FUTURE WORK 


• Assume a BC on the upper surface of the material, x = x 0 

• Using the direct problem computational formula (one-dimensional, implicit, 
centered, finite volume method), solve for thermal solution between the material 
surface and lower BC using x 0 and x 2 as BC's 

• Adjust x 0 estimation until thermal solution at x = x ± matches TC measurement 


estimated BC, x = x 0 Material 


J 

upper BC, x = x t 

i. 

Direct 



lower BC, x = x 2 


TC plug 



Lower TC 


Upper TC 
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